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shown that a vorticity-velocity formulation exists which has the interesting 
property that non-inertial effects only enter into the problem through the 
implementation of initial and boundary conditions. This valuable characteris- 
tic, along with other advantages of the vorticity-velocity approach, are 
discussed in detail. 
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Two distinctly different approaches have been utilized in the literature 
for the numerical solution of the equations of viscous flow in three- 
dimensions. In the more common approach, the momentum equation, which 
contains both the velocity and pressure, is solved numerically along with a 
derived Poisson equation for the pressure (i.e., the pressure-velocity or 
primitive variable formulation [1-3]). The alternative approach is based on 
eliminating the pressure from the momentum equation by the application of the 
curl. In this manner, a vorticity transport equation is solved numerically in 
lieu of the momentum equation (i.e., the vorticity-velocity formulation [4- 
6]). The purpose of the present note is to explore in more detail the 
properties of these disparate numerical approaches. It will be shown that the 
vorticity-velocity formulation has a striking advantage when applied to 
problems in non-inertial frames of reference. More specifically, there exists 
an intrinsic vorticity-velocity formulation wherein all non-inertial effects 
(arising from both the rotation and translation of the frame of reference 
relative to an inertial framing) only enter into the solution of the problem 
through the implementation of initial and boundary conditions . This is in 
stark contrast to the pressure - velocity formulation where non-inertial 
effects appear directly in the momentum equation in the form of Coriolis and 
Eulerian accelerations — a state of affairs which can give rise to a variety of 
numerical problems [2]. A detailed exposition of this interesting property of 
the vorticity-velocity formulation will be presented along with a brief 
discussion of other advantages of this approach. 

For simplicity, we will restrict our attention to the analysis of viscous 
incompressible flow governed by the Navier-Stokes equation and continuity 
equation which, respectively, take the form 
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3v 

— + v • Vv = - Vp + vV“ 
d t ~ — ~ 


( 1 ) 


V • v = 0, 


( 2 ) 


where v is the velocity vector, p is the pressure, and v is the kinematic 
viscosity of the fluid. Here, the validity of (1) requires that the external 
body forces be conservative and that the frame of reference be inertial. In 
an arbitrary non-inertial frame of reference (which can rotate with a time- 
dependent angular velocity fi(t) and translate with a time-dependent velocity 
Vg(t) relative to its origin 0), the Navier-Stokes equation takes the more 
complex form [7] 

3v # 2 
- — + v*Vv+ftxr+fix (fl x r) + v. + 21! x v = - Vp + vV v. (3) 

Here, r is the position vector and the non-inertial terms on the left-hand 
side of (3) are, respectively, referred to as the Eulerian, centrifugal, 
translational, and Coriolis accelerations. The continuity equation still 
assumes the same form (2) in any non-inertial frame of reference. 

By the introduction of a modified pressure P which includes the 
centrifugal and translational acceleration potentials, the non-inertial form 
of the Navier-Stokes equation (3) can be simplified considerably. More 
specifically, (3) can be written in the equivalent form 

3v # 2 

+ v*Vv+Bxr+22xv = -VP + vV v, (4) 


where 
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P = p + 


7 (3 ‘ £) 2 - T 


2 r 2 + 


r. 


(5) 


In the pressure-velocity formulation, equation (4) is solved in conjunction 
with a Poisson equation for the pressure which is obtained by taking the 
divergence of (4). Hence, the governing equations to be solved numerically in 
this approach can be summarized as follows: 

3v ^ 2 

- — +v*Vv+ftxr+2fixv=-VP+vV v, (6) 

dt~~~~~ ~ ~ ~ ~ 

V 2 P = - trfVv • Vv) + 2ft • a), (7) 

' fssrxs y r*j /v» 


subject to the initial and boundary conditions 
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A/ 
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0 * 


at t 


t 


0 ’ 


( 8 ) 


v 

P 



on B. 


(9) 


In (7) and (9), tr(*) denotes the trace, m is the vorticity vector, and B 
denotes the boundary surface of the region. Of course, equations (6) and (7) 
must be solved subject to the continuity equation (2). Since we are 
considering general three-dimensional flow, a stream function solution does 
not exist. Hence, the solution for the velocity v must be projected in some 
suitable fashion onto the space of solenoidal vectors. 

It is quite clear that the form of (6) and (7) (and, hence, their 
mathematical character) change depending on whether or not the frame of 



-4- 


reference is inertial. Consequently, a particular numerical algorithm which 
may he optimal for a given class of flows in an inertial frame of reference 
may not be so for the same class of flows in a non-inertial framing. It will 
now be demonstrated that the vorticity-velocity formulation does not suffer 
from this deficiency. 

The vorticity-velocity formulation is based on the vorticity transport 
equation which is obtained by taking the curl of (4). This equation takes the 
form 

3u „ 

t — + v • Vw = to • Vv + vV u + 20 , • Vv - 2ft (10) 

d f 'V /v<v /w«-w aw /ww rw 

in any non-inertial frame of reference where 


m 


V x v 


(ID 


is the vorticity vector. It is clear that the velocity and vorticity are also 
connected through the Poisson equation 


V 2 v 


= - V x 


( 12 ) 


which is a direct consequence of the vector identity 


V x (V x v) = V(V • v) 

<w rs* ~ rw <-w /w 


V. 

/w 


(13) 


The intrinsic vorticity W, defined by 


y 


(14) 


W = a) + 20 . 
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can be introduced which represents the vorticity relative to an inertial frame 
of reference. Since is spatially homogeneous (i.e., V£2 = 0), it is a 

simple matter to show that the non-inertial form of the vorticity-velocity 
formulation can be written as follows: 

3W 

_21 + v • VW = W • Vv + vV W (15) 

V 2 v = - V x W. (16) 


Equations (15) - (16) must be solved (in some region R with a boundary 
surface B) subject to the initial and boundary conditions 


W = (V x v) Q + 2£ 0 , 


at t 


(17) 


v = v T 


on B. 


W = (V x v) + 20 
~ ~ ~ B ~ 


(18) 


Of course, it is well known that the vorticity, as well as the intrinsic 
vorticity, are solenoidal, i.e., 


V • W = 0, (19) 

~ rst 

and, hence, the solutions for W and v must, in some suitable fashion, be 
projected onto the space of solenoidal vectors. 

This vorticity-velocity formulation of fluid dynamics represented by 
equations (15) - (18) has the striking property that non-inertial effects only 
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enter into the solution of the problem through the implementation of initial 
and boundary conditions . Consequently, the basic structure of the numerical 
algorithm (i.e., the numerical formulation of (15) - (16)) will be independent 
of whether or not the frame of reference is inertial — a situation which 
greatly enhances the general applicability of any Navier-Stokes computer code 
which is developed based on this approach. 

At this point , a few comments should be made concerning the alternate 
ways in which the velocity field can be calculated in the vorticity-velocity 
formulation. Instead of solving the Poisson equation (16), it is possible to 
solve the defining equation for vorticity directly, i.e., 

Vxv=m=W- 2 2, (20) 

(see Gatski, Grosch, and Rose [6]). Of course, for plane or axisymmetric 
flows, there exists a stream function ip such that [7] 

v = A x Vil; (21) 

/v iv 1 

V x (A x V*) = W - 20, (22) 

where A = Vy and y is the coordinate that the flow is independent of (for 
plane flows, (22) reduces to the Poisson equation V = W - 20). While the 

motion of the frame of reference does enter into the equations of motion in 
these alternate vorticity-velocity formulations, it does so in a much less 
significant way than in the pressure-velocity formulation. To be specific, 
the transport equation which is solved (i.e., equation (15)) does not contain 
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any frame-dependent terms and, at each time step, the partial differential 
equation for the determination of the velocity field is only altered by the 
addition of a constant forcing function in the form of 2£2 (the added term on 
the right-hand side of (20) and (22)). 

Finally, it would be of value to mention some other advantages of the 
vorticity-velocity formulation. More difficulties have been known to arise in 
the implementation of pressure boundary conditions than vorticity boundary 
conditions [1,2] (of course, both boundary conditions must usually be 
derived). Difficulties in satisfying the continuity equation in the pressure- 
velocity formulation have also been known to give rise to numerical 
instabilities [1], Furthermore, in the vorticity-velocity approach, the 
vorticity vector is calculated directly. This is of considerable value since 
the vorticity field can play an important role in characterizing certain 
features of turbulence [8], While it is certainly not being suggested that 
the pressure-velocity formulation be abandoned, this study does indicate that 
the vorticity-velocity formulation can have distinct advantages when applied 
to an important class of viscous flows. 
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